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Abstract 

The Harary index of a graph G is recently introduced topo- 
logical index, defined on the reverse distance matrix as H{G) — 
vev(G) d{u v) ' where d{u, v) is the length of the shortest path be- 
tween two distinct vertices u and v. We present the partial ordering 
of starlike trees based on the Harary index and we describe the trees 
with the second maximal and the second minimal Harary index. In 
this paper, we investigate the Harary index of trees with k pendent 
vertices and determine the extremal trees with maximal Harary in- 
dex. We also characterize the extremal trees with maximal Harary 
index with respect to the number of vertices of degree two, matching 
number, independence number, radius and diameter. In addition, 
we characterize the extremal trees with minimal Harary index and 
given maximum degree. We concluded that in all presented classes, 
the trees with maximal Harary index are exactly those trees with the 
minimal Wiener index, and vice versa. 
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1 Introduction 

In theoretical chemistry molecular structure descriptors (also called topo- 
logical indices) are used for modeling physico-chemical, pharmacologic, tox- 
icologic, biological and other properties of chemical compounds [111 [24] , 
There exist several types of such indices, especially those based on graph- 
theoretical distances. 

In 1993 Plavsic et al. in j22_ and Ivanciuc et al. in [T5] independently in- 
troduced a new topological index, which was named Harary index in honor 
of Frank Harary on the occasion of his 70th birthday. This topological 
index is derived from the reciprocal distance matrix and has a number of 
interesting chemical-physics properties [TB]. The Harary index and its re- 
lated molecular descriptors have shown some success in structure-property 
correlations [H [31 SI [TUl HZl HOI HH] ■ Its modification has also been proposed 
|21| and their use in combination with other molecular descriptors improves 
the correlations '53] . It is of interest to study spectra and polynomials 
of these matrices [9, 27^. 

In this paper, let G be a simple connected (molecular) graph with vertex 
set V{G). The Harary index is defined as the half-sum of the elements in 
the reciprocal distance matrix (also called the Harary matrix [18]), 



where d(u, v) is the distance between u and w in G and the sum goes over 
all the pairs of vertices. The Wiener index, defined as 



is considered as one of the most used topological indices with high correla- 
tion with many physical and chemical properties of molecular compounds. 
The majority of chemical applications of the Wiener index deal with acyclic 
organic molecules. For recent results and applications of Wiener index 
see [S]. 

Up to now, many results were obtained concerning the Harary index of 
a graph. Gutman [101 supported the use of Harary index as a measure of 
branching in alkanes, by showing 

Theorem 1.1 Let T be a tree on n vertices. Then 



The right equality holds if and only if T = Sn, while the left equality holds 
if and only if T P„ ■ 




W{G)= '^("'^)' 



u,v£V{G) 




k=2 
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In this paper, we further refine this relation by introducing long chain 
of inequalities and obtain the trees with the second maximum and the 
second minimum Harary index. In |29) Zhou, Cai and Trinajstic presented 
some lower and upper bounds for the Harary index of connected graphs, 
triangle-free and quadrangle-free graphs and gave the Nordhaus-Gaddum- 
type inequalities. In [T] the authors obtained some lower and upper bounds 
for the Harary index of graphs in terms of the diameter and the number of 
connected components. Zhou, Du and Trinajstic j30| discussed the Harary 
index of landscape graphs, which have found applications in ecology [?]. 
Feng and Ilic in [7] establish sharp upper bounds for the Zagreb indices, 
sharp upper bound for the Harary index and sharp lower bound for the 
hyper- Wiener index of graphs with a given matching number. 

In this paper, we analyze relations between the extremal trees with 
maximum and minimum Harary and Wiener index. We investigate the 
Harary index of n- vertex trees with given number of pendent vertices and 
determine the extremal trees with maximal H{G). Furthermore, we derive 
the partial ordering of starlike trees based on the majorization inequalities 
of the pendent path lengths. We characterize the extremal trees with max- 
imal Harary index in terms of the number of vertices of degree two, the 
matching number, independence number, radius and diameter. Finally, we 
characterize the extremal trees with minimal Harary index with respect to 
the maximum vertex degree. All these results are compared with those of 
the ordinary Wiener index. We conclude the paper by posing a conjecture 
regarding to the extremal tree with the maximum Harary index among 
n-vertex trees with fixed maximum degree. 

2 Preliminaries 

For any two vertices u and v in G, the distance between u and v, denoted 
by dciu, v), is the number of edges in a shortest path joining u and v. The 
eccentricity e{v) of a vertex v is the maximum distance from v to any other 
vertex. The vertices of minimum eccentricity form the center. A tree has 
exactly one or two adjacent center vertices; in this latter case one speaks of a 
bicenter. The diameter d{G) of a graph G is the maximum eccentricity over 
all vertices in a graph, while the radius r{G) is the minimum eccentricity 
over all v G V{G). For a vertex u in G, the degree of u is denoted by 
deg{u). 

Two distinct edges in a graph G are independent if they are not incident 
with a common vertex in G. A set of pairwise independent edges in G is 
called a matching in G, while a matching of maximum cardinality is a max- 
imum matching in G. The matching number P[G) of G is the cardinality of 
a maximum matching of G. It is well known that /3(G) < with equality 
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if and only if G has a perfect matching. The independence number of G, 
denoted by a{G), is the size of a maximum independent set of G. 

Let Pn and 5„ denote the path and the star on n vertices. A starlike 
tree is a tree with exactly one vertex of degree at least 3. We denote by 
S{ni,n2, . . . ,nk) the starlike tree of order n having a branching vertex v 
and 

S{ni,n2,...,nk) - = P„i U P„2 U . . . U P„^, 

where ni > n2 > . . . > Uk > 1. Clearly, the numbers ni, n2, . . . , «.fe 
determine the star like tree up to isomorphism and n = ni+n2 + ■ ■ ■+nk + l- 
The starlike tree BSn,k — S{n\,n2, ■ ■ ■ ,nk) is balanced if all paths have 
almost equal lengths, i.e., jn, — nj| ^ 1 for every 1 S^i < j ^ k. 

Denote by A(r) the maximum vertex degree of a tree T. The path P„ 
is the unique tree with A = 2; while the star Sn is the unique tree with 
A = n — 1. Therefore, we can assume that 3 < A < n — 2. The broom 
Bn,A is a tree consisting of a star S'a+i and a path of length n — A — 1 
attached to an arbitrary pendent vertex of the star. 

If ^^^Y^ < m < n— 1, then An,m is the tree obtained from Sm+i by adding 
a pendent edge to each of n — m — 1 of the the pendent vertices of Sm+i- We 
call An,m a spur (see Fig. 1). Clearly, An,m has n vertices and m pendent 
vertices; the matching number and the independence number of An,m are 
n — m and m, respectively. Note that if to > then An^m — BSn,m- 




Figure 1: The spur ^13,7 . 

Next, we give some lemmas which are very useful in the following. 

Definition 2.1 Let v be a vertex of a tree T and deg{v) — m + 1. Sup- 
pose that Pi, P2, • • • , Pm are pendent paths incident with v, with the start- 
ing vertices of paths vi , ?;2 , . . . , Vm , respectively and lengths nt ^ 1 (i = 
1,2,..., to). Let w be the neighbor of v distinct from Vi. Let T' = 5{T, v) 
be a tree obtained from T by removing the edges vvi,vv2, ■ ■ ■ ,vvm-i and 
adding edges wvi , WV2, ■ ■ ■, wvm-i- We say that T' is a S -transform ofT. 

This transformation preserves the number of pendent vertices in a tree T. 
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Figure 2: 5 -transformation on vertex v. 

Lemma 2.2 Let T be a tree rooted at the center vertex u with at least 
two vertices of degree 3. Let v G {z\ deg{z) > 3,z ^ u} be a vertex 
with the largest distance d{u, v) from the center vertex. Then for the S- 
transformation tree T' — S{T,v), it holds 

H{T') > H{T). 



Proof. We follow the symbols in Definition 12.11 Let G be the com- 
ponent of T — wv containing the vertex w (as shown in Fig. 3). Let 
R = {Pi, P2, . . . , Pm-i}- After (5-transformation, the distances between 
vertices from G and R decreased by 1, while the distances between vertices 
from R and Q = Pm U {v} increased by 1. By direct calculation, we have 



H{T')-H{T) = ^ 



1 



1 



d(x,y) — l d(x,y) 



E 



1 



1 



xeQ.yeR 



d(x,y) + l d{x,y)' 



According to the assumption, there is an induced path P ~ wwiW2 



■Wk 



in G, with length at least max{ni, ^2, . . . , rim}. For each path Pi, 1 < i < 



m 



1, it follows 



D, 



xGG,yePi 

E 

x&P,yePi 

> 0. 



1 



> 



d{x,y) 
1 



1 



d{x,y) 
1 



d{x,y)-l d{x,y) 



1 



1 



^ \d{x,y) + l d(x,y) 
xeQ,yeP, ^ ^ V 



1 



1 



f-^ „ \d(x,y) d{x,y) + 1 
x£Q,y£Pi \ V V ly; I 
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Hence 

m — l 

H{T') - H{T) = ^ A > 0. 

1=1 

Since T contains at least two vertices of degree at least 3, we have strict 
inequality. Therefore, if we move pendent paths Pi towards the center 
vertex m of T along the path P, the Harary index increases. I 

Lemma 2.3 U0\/ Let G be a connected graph and v £ V{G). Suppose 
that P = vqVi . . . Vk, Q — uqUi . . . Um cire two paths with lengths k and m 
(k > m > 1), respectively. Let Gk,m. be the graph obtained from G, P, Q by 
identifying v with vq, u with uq, respectively. Then 

H{Gk,m) > H{Gk+l,m-l)- 

By Lemma 12.21 and Lemma 12.31 we can get the following 

Proposition 2.4 Let Gq be a connected graph and u £ V{Go). Assume 
that Gi is the graph obtained from Go by attaching a tree T (T ^ Pk and 
T ^ Sk) of order k to u; G2 is the graph obtained from Gq by identifying 
u with an endvertex of a path Pk; G3 is the graph obtained from Gq by 
identifying u with the center of a star Sk. Then 

H{G2)<H{Gi) < HiGs). 

Let X = {xi,X2, . . . , x„) and y = (yi, y2, • ■ • , Vn) be two integer arrays 
of length n. We say that x majorizes y and write a; )~ y if the elements of 
these arrays satisfy following conditions: 

(i) xi ^ X2 ^ . ■ . ^ Xn and 2/1 > j/2 ■ • • ^ ?/«, 

(ii) xi + X2 + . . . + Xk yi + y2 + ■ ■ ■ + yk, ioT every 1 ^ k < n, 

(iii) xi + X2 + . . . + Xn ^ yi + y2 + ■ ■ ■ + yn- 

Theorem 2.5 Let p — (pi,p2, • ■ • ,Pk) cind q = {qi, qi, . . . , qk) be two ar- 
rays of length k ^ 2, such that p < q and n = pi + p2 + . . . + Pk ~ 
qi + q2 + ■ ■ - Qk- Then 

H{S{pi,P2,...,Pk))>H{Siqi,q2,...,qk)). (1) 

Proof. We will proceed by induction on the size of the array k. For k = 2, 
we can directly apply transformation from Lemma 12.31 on tree 5(91,^2) 
several times, in order to get S{pi,p2). Assume that the inequality ^ holds 
for all lengths less than or equal to fc. If there exist an index 1 ^ m < fc 
such that P1+P2 + . . . +Prn = 91 + 92 + • ■ • + ^Zm, we can apply the induction 
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hypothesis on two parts S{qi,q2,..., q^) U S{qm+i,q,n+2, • • ■ , 9fc) and get 
S{pi,p2, ■ . ■ ,p„i) U S{p 

Otherwise, we have strict inequahties pi+p2 + - • ■+Pm < 'Zi+'72 + - ■ ■+qm 
for all indices 1 ^ m < fc. We can transform tree S{qi,q2, ■ ■ ■ ,qk) into 

S{qi,q2, . . .,qs-i,qs - ^-.qs+l, ■ ■ • ,q,-i,(7r + l,<7r+l, • ■ ■,qk), 

where s is the largest index such that qi = (72 = • ■ • = <?s and r is the 
smallest index such that qr = qr+i = ■ ■ ■ = qk- The condition p ~< q is 
preserved, and we can continue until the array q transforms into p, while 
at every step we increase the Harary index. I 

Corollary 2.6 Let T = S{ni, n2, . . . , n-k) be a starlike tree with n vertices 
and k pendent paths. Then 

HiBn,k) < H{T) < H{BSn,k)- 

The left equality holds if and only if T = Bn.k o,nd the right equality holds 
if and only if T BSn,k- 

3 Main results 

3.1 Trees with given number of pendent vertices 

Let %i,k (2 < fc < n — 1) be the set of trees on n vertices with k pendent 
vertices. If fc = 2, then the tree is just the path P„. If fc = n — 1, then 
the tree is just the star Sn- Therefore, we can assume that 3 < fc < n — 2 
in the sequel. It was proved in [13] that among n-vertex trees with given 
number fc of pendent vertices or given number q of vertices of degree two, 
BSn.k and BSn.n-q-q havc minimal Wiener index. 

Theorem 3.1 Of all the trees on n vertices with k {3 < k < n — 2) pendent 
vertices, BSn.k is the unique tree having maximal Harary index. 

Proof. Suppose T G %i,k has the maximal Harary index, rooted at the 
center vertex. Let St — {v ^ y{T) ■ deg{v) > 3}. 

If \St\ = 1, then by Corollarv l2.6l it follows that BS{n, fc) is the unique 
tree that maximizes Harary index. If \St\ > 2, then there must exist at 
least two vertices of degree at least 3 and there are only pendent paths 
attached below them. We can consider T as the rooted tree at the center 
vertex, and choose the vertex v of degree at least 3 that is furthest from 
the center vertex. After applying (5-transformation, we increase H(G) while 
keeping the number of pendent vertices fixed - which is contradiction. I 
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Theorem 3.2 Among trees with fixed number 9 (0 < g < rt — 1) of vertices 
of degree two, Tn^n-i-q is the unique tree having maximal Harary index. 

Proof. The proof is similar to that of Theorem 13.11 

Consider (5-transformed tree T' — S{T,v). The vertex v has degree 
greater than two, while the vertex w has degree greater than or equal to 
two. Among vertices on the pendent paths Pi, P2, ■ ■ ■ , Pm, there are 

m 

S" = (ni - 1) + (n2 - 1) + . . . + (n,„ - 1) = ^ - m 

i=l 

vertices of degree two. 

If w has degree two in T, then after (5-transformation, v will have degree 
two, and the number of vertices of degree two in T' remains the same as in 
T. Otherwise, assume that deg{w) > 2. We can apply one transformation 
from Lemma 12. 3[ and get new tree T" with m + 1 pendent paths attached 
at vertex w with lengths ni, 712, . . . , Um, 1- This way we increased Harary 
index, while the number of vertices with degree two in trees T and T" are 
the same. 

By repetitive application of these transformations, we prove that the 
starlike tree has maximal value of Harary index among trees with q vertices 
of degree two. The number of pendent paths is exactly k = '^j ~ 

q = n — 1 — q, and by CoroUarv 12.61 it follows that balanced starlike tree 
BS{n, n — 1 — q) is the unique tree that maximizes Harary index. I 

By Lemma 12.31 we have the following chain of inequalities 

H{P„) = H{BSn,2) < H{BSn,^) <■■■< H{BSn,n-i) = H{Sn). (2) 

Notice that Lemma 2.2 from Hj, iJ(F„) < ("+2K"-i) ^ H{Sn) follows 
directly from these inequalities. 

3.2 Trees with given matching or independence num- 
ber 

Du and Zhou in ^ proved that among n-vertex trees with given matching 
number m, the spur An,m minimizes the Wiener index. 

Lemma 3.3 For arbitrary 2^ < m < n — 1, there holds 

H{An,m) = ^ (3^2 + 2mn + - 9m + 19n - 22) . 

Proof. There are four types of vertices in the tree An^m- Denote with 
D'iv) the sum of all reverse distances from v to all other vertices. 
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m I n—m — 1 . 



For the center vertex, D'{v) — ^ 

For each pendant vertex attached to the center vertex, D'{v) = j 



2 



m— 1 _|_ n—m — 1. 



• For each vertex of degree 2, different from the center vertex, D'{v) = 

2 I m— 1 I n—m — 2 . 
1 2 3 ' 

• For each pendant vertex not attached to the center vertex, D'{v) = 

1 ' 2 ^ 3 ' 4 

After summing above contributions to the Harary index, we get 

n + m — I , ,2ri + m + l 

i?(An.m) = 1 + (2m-n+l) 

z 

2ri + m + 5 , ,3n + m + 8 
+ (n - m - 1) • g h (n - m — 1) • 

= — iir? + -22- 9m + I9n + 2mn) . 
24 ^ ^ 

This proves the resuhs. I 



Theorem 3.4 Let T be a tree on n vertices with matching number /?. Then 

H{T) < ^ {6n^ - 4:(3n + /3^ +9(3 + lOn - 22) , 

with equality holding if and only if T ^ An.n-p- 
Proof. Suppose that T has k pendent vertices, then 

fc</3 + n-2/3 = n-/?, 

and by Theorem 13. 1[ we have H{T) < H{BSn,k)- Using equation (2)), it 
follows H{BSn,k) < H{BSn,n-p) = 7?(A„,„_^)' since n - /3 > f . Finally, 
H{T) < H{An^n-i3), with equality holding if and only if T = An^n-p- By 
Lemma liSl we obtain the explicit relation for H (^An^rt—p^ • B 

By Theorem 13.41 we have the following corollary. 
Corollary 3.5 Let T be a tree of order n with perfect matching. Then 

H{T) < i(17n2 + 58n-88), 
with equality holding if and only ifT'^ An,zL. 
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Theorem 3.6 Let T be a tree on n vertices with independence number a. 
Then 

H{T) < — {3n^ + 2an + ~ 9a + 19n - 22) , 
with equality holding if and only if T ^n.a- 

Proof. Since all pendent vertices form an independent set, it follows k < a. 
Every tree is bipartite graph, and we get a > \^~\ . By Theorem l3.11 we have 
H{T) < H{BSn,k) and H{BSn,k) < H{BSn,c) = i?(A„,„). Therefore, 
H{T) < H{An^a), with equality holding if and only if T = An^a- H 

3.3 Trees with given diameter or radius 

Let C'n,d{pi,P2, ■ ■ ■ ,Pd-i) be a caterpillar on n vertices obtained from a 
path Pd+i — vqVi . . -Vd-iVd by attaching Pi > pendant vertices to w^, 
1 < i < d — 1, where n = d + 1 + X^iLi^ Pi- Denote 

Cn,d,i = Cn,d{0, . . . , 0, n - - 1, . . . , 0). 

Obviously, Cn,d.i — Cn,d.n-i- In [14| and |19| it is shown that caterpillar 
Cn,d.ld/2\ bas maximal Wiener index among trees with fixed diameter d. 

Theorem 3.7 Among trees on n vertices and diameter d, C'n,d,[d/2i is the 
unique tree having maximal Harary index. 

Proof. Let T be n-vertex tree with diameter d having maximal Harary 
index. Let Pd+i = vqVi ... be a path of length d. By Proposition 12.41 
all trees attached to the path Pd+i must be stars, which implies that T = 
Cn,d{Pi,P2, ■ ■ ■ ,Pd-i)- Applying Lemma \T2\ at those vertices of degree at 
least 3 in Cn,d{pi,P2, ■ ■ ■ ,Pd-i), we get that T = Cn,d,i for some 1 < i < 
d — 1. By Lemma 12.31 we get that the maximal Harary index is achieved 
uniquely for i = [^J. This completes the proof. I 

For a tree T with radius r(r), it holds 2r{T) = diT) or 2r(T)-l = d{T). 
Using transformation from Lemma |2 . 31 applied to a center vertex, it follows 

that H{Cn.2r.,yd/2\) < H{Cn,2r-l.ld/2\)- 

Corollary 3.8 Let T be a tree on n vertices with radius r > 2. Then 

ff(T) <ff(C„,2r-i,Ld/2j), 
with equality if and only if T ^ C'n.2r-i.[(i/2j ■ 
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If d > 2, we can apply the transformation from Lemma [^751 at the center 
vertex in C„,rf^Ld/2j to obtain C„^d-i,L(d-i)/2j • Thus, 

H{Pn) ~ H{Cn^n-l,y(rL-l)/2\) < •■• < H{Cn3,l) < H{Cn,2s) — H{Sn)- 

Also, it follows that H{Cn,3.i) has the second maximum Harary index 
among trees on n vertices. 

3.4 Trees with given maximum vertex degree 

Chemical trees (trees with maximum vertex degree at most four) provide 
the graph representations of alkanes jll] . It is therefore a natural problem 
to study trees with bounded maximum degree. It was proven in [5 that 
among n-vertex trees with the maximum degree A, the broom i3„_A has 
maximal Wiener index. 

Lemma 3.9 For arbitrary 2 < A < n — 1, there holds 

^ , , , (A-l)(A-2) , A-1 



- A + 1' 



H{B„^a) = • ffn-A - n + . 

'i II 

where Hk = 1 + ^ + ... + ^ he the k-th harmonic number. 

Proof. The Harary index of i?„,A can be calculated as the sum of Harary 
index of P„_a+i, the sum of reverse distances between A — 1 pendent 
vertices and vertices from a long path, and the sum of reverse distance 
between pendent vertices. 

/A-n 

i?(S„,A) - i/(P„_A+l) + + (A - l)i?„_A + l 



= n ■ H„/ 



(A-l)(A-2) A-1 



A + 1 



Theorem 3.10 Let T be a tree on n vertices with the maximum degree A. 
Then H[T) < H{Bn,A)- The equality holds if and only if T = Bn,A- 

Proof. Fix a vertex v of degree A as a root and let Ti , T2 , . . . , Ta be the 
trees attached at v. By Proposition 12.41 all subtrees attached to Ui are 
paths for 1 < i < A, the Harary index increases. This implies the result. 
■ 

If A > 2, we can apply the transformation from Lemma 12.31 at the 
vertex of degree A in i?„^A and obtain -B„,a-i- Thus, 

HiSn) = HiBn.n-l) > H{Bn,n-2) > ■■■> H{Bn,3) > H{Bn,2) = H {Pn) . 

Also, it follows that i3„,3 has the second minimum Harary index among 
trees on n vertices. 
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4 Concluding remarks 



In this paper, we presented the partial ordering of starhke trees based 
on the Harary index and we derived the trees with the second maximal 
and the second minimal Harary index. We characterized the extremal 
trees with maximal Harary index and fixed number of pendent vertices, the 
number of vertices of degree two, matching number, independence number, 
radius and diameter. In addition, we characterized the extremal trees with 
minimal Harary index and given maximum degree. We concluded that in 
the all presented classes, the trees with maximum values of Harary index 
are exactly those trees with the minimal Wiener index W{G), and vice 
versa. 

The complete A-ary tree is defined as follows. Start with the root having 
A children. Every vertex different from the root, which is not in one of 
the last two levels, has exactly A — 1 children. In the last level, while not 
all nodes have to exist, the nodes that do exist fill the level consecutively. 
Thus, at most one vertex on the level second to last has its degree different 
from A and 1. 

In [12 the authors proposed these trees to be called Volkmann trees, 
as they represent alkanes with minimal Wiener index [8]. The computer 
search among trees with up to 24 vertices reveals that the complete A- 
ary trees attain the maximum values of H{G) among the trees with the 
maximum vertex degree A. 

Conjecture 4.1 For any k ^ 2, the complete A-ary tree has maximum 
value of H{G) among trees on n vertices with maximum degree A. 

It would be interesting for further research to consider the extremal 
unicyclic and bicyclic graphs with respect to Harary index, and compare 
the results with those for the Wiener index. 
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